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Abstract 

In this study we report a physical analysis of the membrane mechanics affecting the size of the highly curved 
region of a lipid nanotube (LNT) that is either connected between a lipid bilayer vesicle and the tip of a glass 
microinjection pipette (tube-only) or between a lipid bilayer vesicle and a vesicle that is attached to the tip of a 
glass microinjection pipette (two-vesicle). For the tube-only configuration (TOC), a micropipette is used to pull a 
LNT into the interior of a surface-immobilized vesicle, where the length of the tube L is determined by the 
distance of the micropipette to the vesicle wall. For the two-vesicle configuration (TVC), a small vesicle is inflated at 
the tip of the micropipette tip and the length of the tube L is in this case determined by the distance between 
the two interconnected vesicles. An electrochemical method monitoring diffusion of electroactive molecules 
through the nanotube has been used to determine the radius of the nanotube R as a function of nanotube length 
L for the two configurations. The data show that the LNT connected in the TVC constricts to a smaller radius in 
comparison to the tube-only mode and that tube radius shrinks at shorter tube lengths. To explain these 
electrochemical data, we developed a theoretical model taking into account the free energy of the membrane 
regions of the vesicles, the LNT and the high curvature junctions. In particular, this model allows us to estimate the 
surface tension coefficients from R{L) measurements. 



Background 

Membrane tethers have been studied extensively over 
the past 40 years [1-11]. These structures, also called 
membrane nanotubes, were observed during fluid shear 
deformation of live cells attached to a substrate. As 
these cells were dislodged, membranous tethers 
remained attached to the surface displaying both the 
fluid and the elastic properties of the membrane [1,2]. 
Following this work many naturally forming membrane 
nanotubes have been identified [7-10]. For example, 
membrane nanotubes have been shown to exist within 
the cell, notably in the trans golgi network [10]. Here, 
lipid and protein cargo destined for various destinations 
throughout the cell are sorted and pinched off from the 
tubular membrane of the network. It has also been 
reported that cells have the ability to use membrane 
nanotubes for the exchange of organelles [7], and this 
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exchange has interestingly even been recognized 
between different cell types [8]. Thus, these tethers, 
which were first observed following a dramatic manipu- 
lation, have been shown to be a common occurrence in 
biology. 

Following their initial discovery, the lipid membrane 
nanotubes (LNTs) have been created artificially in sev- 
eral model membrane systems. By attaching a bead or a 
micropipette to a point on the membrane and applying 
a localized mechanical force to the bilayer surface it has 
been shown that a lipid tether can be pulled from the 
vesicle membrane [3-5,11]. The size of the structure is a 
result of the interplay between the curvature elasticity 
effects maintaining the original geometry and the mem- 
brane tension [12]. Tether pulling experiments can be 
used for estimations of tube diameters. By measuring 
the forces required for pulling a tube, the diameter of 
the LNTs were estimated to be 50-200 nm [13]. From a 
tube coalescence method [14] and video pixel analysis of 
accumulated fluorescence images as well as from micro- 
graphs obtained with differential interference contrast 
optics [5], the LNT diameters were determined to be in 
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the range of 100-300 nm [13]. To complement these 
methods, we developed an electrochemical method to 
monitor the diffusion of electroactive molecules through 
the LNT, thus allowing the LNT diameter to be mea- 
sured as a function of nanotube length [11]. The 
method relies on the formation of a vesicle-LNT net- 
work by using a micropipette technique [5,15]. The 
micropipette-assisted vesicle-LNT network formation 
allows us to create complex systems of vesicles intercon- 
nected by LNTs, including a so-called inward configura- 
tion where a small daughter vesicle is created inside a 
larger mother vesicle, the two vesicles being connected 
by a LNT [6] (see Figure 1A). During network forma- 
tion, the LNT is pulled with a micropipette to the inter- 
ior of the vesicle and thus the opening of the tube faces 
outward to the exterior of the vesicle. This makes it 
possible to monitor the diffusion of a marker molecule 
from the micropipette, through the tube, and out of the 
nanotube opening. The concentration of the molecules 
measured at the opening of the LNT is directly related 
to the inner diameter of a LNT of determined length 
[11]. In this article we use the electrochemical method 
for monitoring the size of a nanotube attached directly 
to the micropipette in the configuration we refer to as 
the tube-only configuration (TOC) (see Figure IB). 

Additionally, by inflating a small ("daughter") vesicle 
at the tip of the micropipette, the diameter of a nano- 
tube placed in between the inner vesicle and mem- 
brane of the outer vesicle can be examined in a 
configuration here called two-vesicle configuration 
(TVC) (see Figure 1A). The measurements show that 



c ^ 



A Two- Vesicles-Configuration (TVC) 




Figure 1 Experimental configurations. Sketches of the 
geometries of the large unilamellar vesicles interconnected with a 
common LNT; (A) the "two-vesicle" configuration, where the LNT is 
connected between the mother vesicle and a small daughter 
vesicle inside of the mother vesicle, (B) the "tube-only" 
configuration where the LNT is connected between the tip of a 
glass pipette and the giant unilamellar vesicle. 



there is a reduction in tube diameter at shorter length, 
and the effect appears to be more pronounced in the 
TVC. In this work we suggest a geometrical model 
based on direct minimization of the Helfrich's func- 
tional for the system of lipid vesicles linked to a LNT 
via junctions of specific geometry. This new model 
presents a unified quantitative analysis of TOC and 
TVC and explains why the length of the LNT in the 
TVC is twice as high as in TOC for a given radius. 
Furthermore, the model has just two parameters, 
which can be chosen to fit the experimental data on 
monitoring of the size of the LNT. This allows for 
identifying the contribution of the surface tension to 
the free elastic energy of the system. This low-tension 
term has been neglected in the related publication 
[11], where a phenomenological description of the sys- 
tem was suggested and only a qualitative consistency 
with experimental data was obtained. 

Experimentals 

Materials and methods 

Surface-immobilized giant unilamellar soybean lipo- 
somes (SBL) were made from soybean polar lipid extract 
(Avanti Polar Lipids, Alabaster, AL), as previously 
described [5,6,11,15]. An injection pipette pulled with a 
commercial pipette puller (Model PE-21, Narishige Inc., 
London, UK) and was back-filled with a 50 mM catechol 
solution. The pipette was then electro-inserted into the 
unilamellar liposome with the aid of a voltage pulse gen- 
erated relative to a 5 /urn counter electrode (ProCFE 
from Dag an Corp, Minneapolis, MN), which was placed 
on the opposite side of the liposome from the injection 
pipette. Carbon fiber working electrodes were fabricated 
in house and have been described elsewhere [11]. Work- 
ing electrodes were held at +800 mV versus a silver/sil- 
ver chloride reference electrode (Scanbur, Sweden). All 
measurements were made using an Axon 200B potentio- 
stat (Molecular Devices, Sunnyvale, CA). 

Nanotube radius measurements and calculations 

The flux of catechol through the nanotube was mea- 
sured using carbon fiber amperometry. A 5 ftm carbon 
fiber microelectrode was placed at the nanotube-lipo- 
some junction. The nanotube was then either length- 
ened or shortened by manipulating the injection pipette. 
After the new length was obtained, the current was 
allowed to stabilize and was then recorded. This process 
was repeated several times for each liposome resulting 
in a series of electrochemical measurements for tubes of 
different lengths. The electrode was then removed from 
the nanotube-liposome junction and allowed to reach a 
steady current to establish a baseline. The difference in 
measured current for a nanotube versus this background 
together with the length of the nanotube was then used 
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to compute the diameter of the nanotube based on the 
previously derived relationship 



R 



N 



Ai 



nFnD- 



AC 



(1) 



where R is the radius of the nanotube of a given 
length L, Ai is the change in measured current with 
respect to the background, n is the number of moles of 
electrons transferred per mole of redox species (for cate- 
chol, this is equal to 2), F ~ 96 485.34 C/mol is Fara- 
day's constant, D = 7.0 x 10~ 6 cm 2 /s is the diffusion 
coefficient of the selected redox species (catechol). AC is 
the change in concentration of catechol over the nano- 
tube length and is equal to the concentration of electro- 
active species in the pipette assuming that the 
concentration at the electrode surface is zero (in our 
experiments AC = 50 mM). 

The results for the tube radius deduced from the 
simultaneous measurement of electrochemical current 
and the tube length by using formula (1) are presented 
in this study. In comparison with our previous publica- 
tion [11], a wider range of the length L of the tube is 
presented for the TOC configuration. 

Theoretical approach 

The system under consideration 

In the first system (Figure 1A), a mother vesicle contains 
a small daughter vesicle on the inside with a common 
LNT connecting the two compartments. In the second 
case (Figure IB), the lipid tube is pulled to the inside of 
the vesicle and is directly fixed to the tip of the micropip- 
ette. Also, there is a source of lipid attached to the 
mother vesicle wall. The presence of lipid source means 
that the surface tension is low. We model the membrane 
as a two-dimensional surface T. Its free elastic energy 
written in the form of Helfrich functional [16] reads 



-1/ 



(2H - CofdA + a A. 



(2) 



Here H is the mean curvature of the surface, C 0 is the 
spontaneous curvature which is determined by the spe- 
cific chemical composition of the membrane, k is the 
coefficient of membrane bending, a is the coefficient of 
membrane surface tension. The equilibrium shape of 
the membrane with pulled cylindrical tubule can be 
found from minimum of the functional 



(2H- 



CofdA + aA-fL, 



where / is the force needed to pull the lipid tube of 
length L [12]. In the case, when the junctions are not 
taken into account, the interplay between membrane 
bending k and membrane tension a produces variability 
in tubule radius and the force f 0 



R 0 = y/k/2c7,f 0 



(3) 



where f 0 is the force needed to hold the tube of radius 
R 0 at a fixed position [12]. However, it was shown that 
for lipid vesicles interconnected with LNTs, either 
pulled outward from the vesicle wall [5,15] or inward 
into the vesicle interior [11,17], the neck elements (the 
junctions between the lipid tube and the vesicle body) 
also contribute to the total free energy of the mem- 
brane. Below we consider a theoretical model based on 
the Helfrich functional to find the equilibrium shape of 
the membrane accounting for the junctions of the speci- 
fic geometry. By comparing the results of numerical 
computations with experimental data, we are able to 
determine the tension in the LNT after fitting the 
experimental data with the geometrical model described 
below. 

The geometrical model 

When the inner vesicle or the junction between the 
micropipette and the nanotube is subjected to the trans- 
lation movement along the LNT axis, the length of the 
tube is changed (increasing or decreasing its value in a 
controlled way, which can be monitored under the 
microscope). During these manipulations the radius of 
the tube adapts to minimize the Helfrich's free energy 
(2) with C 0 = 0, as we neglect any contribution from 
spontaneous curvature. 

We assume that the shape of LNT can be approxi- 
mated by a cylindrical surface of radius R and length L. 
Since radii of both vesicles are much larger than the 
tube radius, the junctions between the cylinder and vesi- 
cles are modelled by toroidal surfaces with the inner 
radius R + r and crossection radius r (Figure 2). In the 
TOC, when the inner vesicle is not present, only one 
junction is considered. Although the junction between 
the micropipette and the tube contributes to the total 
free energy, it is assumed that this contribution does 
not depend on the tube radius R and, thus, the corre- 
sponding term vanishes after the variation. In these set- 
tings, the radius-dependent part of the free energy is 
given by the expression: 



F(L, R, r,k,cr) 



IS 



(2H) 2 dA + a A 



(4) 
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where f = Vq U Ft an< l 



r c = 



flcos0 I ,x e (O,L),0 e (0,2tt) 
.Rsin<? 

^ r sin 0 

{R + r{l - cos©))cos0 I ,© 6 ((1 - 
K (R + r(l - cos©)) sin 0 / 



) G (0,2;r) 



(5) 



The toroidal part of the surface can be parametrised 
by (5) due to translation invariance of the energy func- 
tional (4). The multiplier v assumes the value 1 for 
TOC and 2 for TVC to represent both junctions. 

In Equation 4, L, r, <7 are fixed parameters while the 
radius of the tube R is adjusted to satisfy 



9F(L, R, r,k,cr) 
3R 



0. 



(6) 



The variation (6) yields the following relation between 
the tube length L and radius R 



L{R)- 



vR 2 



i - 2R 2 a 



(r + R) ({R + rf 



2r + R 



rR 3 / 2 (2r + R) 3/2 



, (7) 



where o = ^"-The model parameter r as well as a are 

chosen to obtain the best fit to the experimental data. 
Assuming that the radius of the tube R is much larger 
than the parameter r, the first two terms of the power 
series expansion for (7) with respect to rlR can also be 
used to quantitatively model the measured relation L(R). 
This, simplified, form of (7) reads 



L(R) 



7tv(2R 2 -r 2 (l-4R 2 &)) 



4r (l - 2R 2 a) 
and allows for expressing R as a function of L 



R(L) 



r (4L + vnr) 



2v7t + Aro (2L + vnr) 



(8) 



(9) 



An important feature of the proposed model is the 
asymptote R 0 = ^k/2cr (compared to (3)), which is pre- 
sent in all three relations (7), (8), and (9). As we 
increase L, the radius R grows and the energy of the 
cylindrical part of the surface becomes dominant over 
the energy of the toroidal junctions. Thus, in the limit 
case L — > «>, we obtain the junction free equilibrium 
value of R given by (3). 

Fitting the parameters 

For given K measurements (L if Rj), i = 1 ... 7<T, we vary d 
and r to minimize 



Gi (a,r)=^(L(^)-L0 2 



(10) 



i=i 



by means of conjugate gradient minimization proce- 
dure. Here, the relation L(R) is given by (7). When the 
ratio rlR is small, the approximation (8) can be used 
instead. In this case, one can also fit (9) to the data by 
minimizing the functional 



G 2 (a,r) = J2( R ( L i)- R if- 



(11) 



The latter method is preferable when the relative mea- 
surement error for R is greater than the one for L. 

Results and discussion 

When fitting the curve (7) to the dataset for the TVC, 
the parameter values are r « 1.7 nm and oik ~ 89 pm 2 . 
The corresponding values for the dataset in the case of 
TOC are r « 1.2 nm and oik ~ 54 ftm 2 . The relation L 
(R) with fitted parameters are plotted on Figure 3 (blue 
curves) together with measured experimental data. As 
expected, the parameter r is much smaller than the 
radius R: rlR <0.06. Therefore, the simplified form (8) 
and its inverse (9) can be used for the given range of 
values of R. Fitting the relation (9) to the measurements 
by minimizing (11) yields oik ~ 98 /urn' 2 , r ~ 1.9 nm 
and oik ~ 72 /urn' 2 , r ~ 1.7 nm for TVC and TOC, 
respectively. The corresponding curves are plotted in 
Figure 3 in red. The model exhibits good agreement 
with the empirical data. A rather large scattering of 
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Figure 3 Comparison of experimental and model results. The measurement points (shown as markers) and the predictions of the model 
(solid lines). Parameters for the model predictions were chosen to minimize functionals (11) (red lines) and (10) (blue lines). 



measurement points at high R values in the TOC case is 
reflected as about 20% difference in parameter values 
when using different approaches to find the best fit. In 
this case, the values obtained through fitting (9), namely 
oik ~ 72 /mi 2 , r ~ 1.7 nm have higher reliability. 

Our model establishes a connection between the data 
from TOC and TVC experiments. It follows directly 
from formula (7), that to reach a given radius R of the 
tube, the length L TVC of the tube in the TVC experi- 
ment must be double of that in TOC arrangement 

L wc (R) = 2L TO c (R) • 

To explore this theoretical prediction, we divide the 
lengths obtained in the experiment with TVC by two and 
plot the resulting data set together with the measure- 
ments for TOC on Figure 4. The optimal parameters of 




Figure 4 Comparison of experimental and model results 
(unified description). The measurement points (shown as markers) 
for both TVC and TOC plotted after dividing the JVC length by two. 
Parameters for the model predictions were chosen to minimize 
functionals (11) (red line) and (10) (blue line). 



the model for this, unified, data are oik ~ 55 ^m" 2 , r ~ 
1.2 nm and elk ~ 71 ftm " 2 , r ~ 1.6 nm for functionals 
(10) and (11), respectively. These values are similar to 
ones for the TOC case since this portion of the data is 
more disperse and has much greater contributions to 
functionals (11) and (10) when compared to the data for 
the TVC case. Figure 4 also shows that the measure- 
ments are in agreement in the region, where they overlap, 
i.e., for values of R between 0.05 and 0.06 ftm. 

Assuming the well established value of bending modu- 
lus k = 10" 12 erg [16], the recalculated coefficients of the 
surface tension are found in the interval of <7 ~ 0.01- 
0.02 dyn/cm. These tension values are much smaller 
comparing to the lipid molecular compressibility (100 
dyn/cm) [18] but much larger than the critical surface 
tension for the instability of the membrane cylinder and 
"pearling" (10~ 5 dyn/cm for (DGDG/DMPC membrane) 
LNTs of radius R ~ 0.3 - 5 ftm found in [19] work) 
while comparable with the magnitude of the lateral ten- 
sion (higher limit) for mutual adhesion of lecithine 
membranes ~10~ 4 erg/cm 2 [20]. 

The small value of the junction radius corresponds to 
the strongly deformed state of the membrane. These 
small values should be considered as order estimates, 
since they are attributes of the assumed toroidal geome- 
try of junctions. The real shape of these junctions is 
probably more complex and, thus, cannot be described 
by just two scalar valued parameters. Although freeze- 
fracture electron microscopy does not reveal bilayers 
with curvature less than 20 nm, the value r ~ 1.5 nm 
which is found from the model is similar to the radius 
curvature of small inverted pores (for example, it is 
known that phospholipids spontaneously form inverted 
membrane structures with the radius varying between 
0.5 and 5 nm, and smallest fusion pores have a calcu- 
lated diameter less than 2.5 nm) [21,22]. 
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Conclusions 

We propose a simple geometrical model for the quanti- 
tative explanation of the experimental results on equili- 
brium geometrical shape and LNTs parameters, R(L), in 
the different configurations. The experimental observa- 
tions show that the nanotube diameter is reduced at 
shorter lengths and also that the diameter is consistently 
smaller for the TVC as compared to the TOC for a 
given length. The observed effect is ascribed to originate 
from the elastic junctions, since the phenomenon is 
accentuated in a system containing two necks connected 
to a vesicle membrane. We approximate the shape of 
these junctions by simple geometrical shapes and 
express the free elastic energy of the membrane in 
terms of the length of the LNT, its radius, the radius of 
the junction and the tension of the membrane. Variation 
of the energy with respect to the nanotube radius yields 
an explicit relation between the radius and the length. 
The relation is in agreement with observed values. The 
model enables estimations of the current surface tension 
coefficient and the curvature at junction regions. The 
estimated values of the surface tension are of order 10 2 
dyn/cm and the curvature value at junctions are com- 
parable to ones at fusion pores. Furthermore, the pro- 
posed model offers a clear explanation of the difference 
in measurements for TVC and TOC: in contrast to 
TOC, the TVC features two junction regions, thus, the 
length of the LNT in this configuration must be twice 
as long to achieve the same value of the radius. 
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